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The Yang-Mills gradient flow in finite volume is used to define a running coupling scheme. As 
our main result the discrete /3 -function, or step scaling function, is calculated for scale change 
s = 3/2 at several lattice spacings for SU (3) gauge theory coupled to Nf — 4 fundamental massless 
fermions. The continuum extrapolation is performed and agreement is found with the continuum 
perturbative results for small renormalized coupling. The case of SU (2) gauge group is briefly 
commented on. 
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1. Introduction 



The Yang-Mills gradient flow is a relatively new addition to the arsenal of non-perturbative 
tools for the study of non-abelian gauge theories [jlj ||, |3|] ; see also [Q, ||] for its use in a slightly 
different context. It has been usefully implemented for a high precision scale determination in 
QCD [||] as well as for a running coupling scheme in a finite volume setup [^], among others. 

All these applications involve the operator TrF^ v F^ v and its derivative with respect to the 
gauge field. The discretized expression of the derivative is then used in a finite step integration 
scheme of the gradient flow, 

dA^t) _ dS YM {A) 
dt dA^ 

The operator TrF^ v F^ v enters the simulations at two other instances. First, it is simply used in the 
gauge action for generating the configurations and second in the observable 
(E(t)) = — 2(TrF^ v (?)F^ v (?) i.e. the field strength tensor squared evaluated at flow time t > 0. 
For the gauge action we use the tree level Symanzik improved action and for E(t) the same clover 
improved discretization is employed as in ^ |(| 0]. For the derivative along the flow ( |1.1| ) we have 
used both the Wilson discretization and the tree level Symanzik variant but found that the scaling 
properties of the latter are more favorable for our present investigation and hence will only use that 
one. 

Our goal is to compute the running coupling in SU (3) gauge theory coupled to Nf = 4 flavors 
of massless fundamental fermions. The main observable is the discrete j5 -function in a finite vol- 
ume setup where the running is with the linear size of the system, similarly to the Schroedinger 
functional ||] and to other finite volume approaches. The discrete j8 -function corresponding to a 
scale change of s is (g 2 (sL) — g 2 (L))/log(s 2 ), where the renormalized coupling g 2 (L) is obtained 
from the gradient flow. 

The Schroedinger functional analysis of the same model, 4-flavor QCD, can be found in [[| 

fill. 



2. Gradient flow scheme 



The gradient flow ( 1 . 1 ) is used to define a renormalized coupling in finite volume by first fixing 
the ratio y/St/L = c and then setting 



mn 2 {t 2 E{t)) 
3(^2 -!)(! + 8(c)) ' 



(2.1) 



where 



8(c) 



C 4 K 2 



+# 4 



-l/c 2 



(2.2) 



is given in terms of the Jacobi elliptic function ~&(q) = E»9 n ■ The numerical factors on the right 



hand side of (2.1) are such that to leading order g 2 -^ = g 2 [0]. Corrections contain both even and 



odd powers of g^. Different choices for the constant c correspond to different schemes but the 



2 



Running coupling 



Daniel Nogradi 



leading order relationship g 2 -^ = g 2 ensures that the 1-loop /3 -function in the gradient flow scheme 
coincides with the 1-loop /3 -function in the MS scheme. The 2-loop /3 -function of the MS scheme 
is on the other hand not the same as in the gradient flow scheme because of the non-zero a\ term: 

£ = 8 m ( 1 + a 1 ( c )#ms + («ms ) ) ; for more details see • 

Massless fermions can also be included if anti-periodic boundary conditions are used. The 

non-trivial boundary conditions cause the fermions to have an effective energy gap of order l/L 

in small volume hence as one follows the running of the coupling from the UV towards the IR at 

least for small renormalized coupling the simulation will not run into problems even at zero bare 

fermion mass. 



3. A note on the SU (2) case 

One might worry that for gauge group SU (2) the matrix integrals in ^ used for the definition 
of the running coupling are not finite in D = 4 dimensions. The matrix integrals needed at leading 
order are 

f ^iTr^fiypexp (iTr[g M ,g v ] 2 ) 

J^exp(lTr[B M ,S v ]2) ' ^ 

where the integral dB is over D anti-hermitian traceless N x N matrices. Throughout this section 
we set L = 1. Both the numerator and denominator are finite in D = 4 and N > 2; see [[TT|]. 
Clearly, one may evaluate Id,n from the matrix model partition function 

Z D , N {b) = J dBexp (\r[B^B v ] 2 ^ (3.2) 

via its logarithmic derivative, 

dlogZ D:N (b) 



1D,N 



db 



(3.3) 

b=l 



provided Zjj.N(b) and its derivative are finite. For D = 4 and N > 2 this is the case and one easily 
obtains I4 jy = N 2 — 1. For D = 4 and N = 2 however Z4^(b) is divergent. But the partition function 
Zo,2{b) can be evaluated for D > 4 and is in fact finite [II! 



7 (f,)- 1 ( 2n2 \ 4 r(D/4)r(D/4-l/2)r(D/4-l) 
D ^ ' 2\ b ) r(D/2)r(D/2-l/2)r(D/2-l) ' K ' 

The divergence in D = 4 is coming from the pole of Y{D/4 — 1). The logarithmic derivative of the 
above expression for D > 4 then leads to 

3D 

Id.2 = -j- , (3-5) 

which agrees with the result Af 2 — 1 for Af = 2 if D = 4 is set at the end of the calculation. Since the 
perturbative calculation of (t 2 E(t)) is performed in dimensional regularization and D = 4 is only 
set at the very end, the above procedure is natural. 
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One needs to be very careful about corrections which are only logarithmically suppressed 



This question was first discussed in Q13|]. A toy model where the issue at hand can be illustrated is 

Z(g 2 ) = / dxdye ' y ^ , (3.6) 
with the associated expectation value {x 2 y 2 }. Let us introduce b = l/g 2 , then 

dZ(b) 

A key feature of the gauge theory case, a flat direction in the tree level potential, is shared by 
the above toy problem. The tree level potential is x 2 y jg and the "1-loop" potential x 2 + y 2 is 
subleading. For finite g both the partition function and the expectation value are finite but the 
g — > limit is quite subtle. 

Naively, in the g — > limit the "1-loop" potential x 2 + y 2 can be dropped to leading order 
but then both the numerator and denominator in are divergent, similarly to the SU (2) gauge 
theory case. 

One might nevertheless first drop the x 2 + y 2 potential and then rescale the variables by x — > 
yjgx and v — > ^fgy in order to obtain 

,22, 2 Jdxdyx 2 y 2 e- x2y2 

{Xy)=8 f dxdye-** + - (3 " 8) 

where again both the numerator and the denominator are divergent. But since now the potential is 
a homogeneous polynomial, the above ratio is naively 1 /2, leading to 

(xY) = £ + ... , (3.9) 

as the naive form of the small-g behavior. However exact evaluation of the integrals for finite g 
leads to 

using the Bessel K-functions. The small-g expansion of the above is 

where y# is Euler's constant. The naively obtained leading order g 2 /2 result is supplemented by 
only logarithmically suppressed 0(g 2 /log(g 2 )) corrections not only polynomial ones like 0(g 4 ). 



Ref. [13] discusses the SU(2) theory, for which a similar logarithmic correction is expected. 
Note that the coefficient of the logarithmic term for the SU(2) theory is at present unknown. 

Since all matrix integrals to leading order for N > 2 are finite we do not expect such complica- 
tions for SU (3) which is our main application. In high orders of g-^ similar logarithmic corrections 
may enter but not to leading order which ensures that in the UV the 1-loop /3 -function is the same 
in our scheme as in every other scheme. 



DN is grateful to Martin Luscher for an instructive discussion of this point. 
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Figure 1: Discrete j3 -function of SU(3) gauge theory coupled to Nf = 4 flavors of massless fundamen- 
tal fermions. Left: scale change of s = 3/2, the results at 3 lattice spacings are shown together with the 
continuum universal 1-loop and MS 2-loop results for comparison. Right: continuum extrapolation of the 
discrete /3 -function for an illustrative value of the coupling g 2 (L) = 3. Two methods were used to perform 
the continuum extrapolation; for details see [^J. The two methods give results which agree with each other 
for all of ourg 2 (L) values. 




Figure 2: Discrete j3 -function of SU (3) gauge theory coupled to Nf = 4 flavors of massless fundamental 
fermions for a scale change of s = 3/2. The continuum extrapolated result is shown together with the 
universal 1-loop and MS 2-loop results for comparison [R]. 



4. Results 



Once the renormalized coupling is defined by (2.1) the discrete j3 -function, or step scaling 
function [[14]], can be computed by increasing the volume by a factor of s. The simulations were 
performed with the same parameters as reported in [j7]|. 

For the scale change s = 3/2 we have results for 4 lattice spacings, corresponding to the 
volume changes 8 — > 12, 12 — > 18, 16 — > 24 and 24 — > 36. The coarse lattice spacing, 8 — > 12, is 
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definitely outside of the a 2 scaling region and can not be used for continuum extrapolation. Thus 
we do not discuss these coarse lattices any further. The results are shown in the left panel of figure 
[I] with the choice of c = 0.3. (Note that as a cross check we repeated the calculations with the scale 
change s = 2 and with 3 lattice spacings, corresponding to 8 — > 16, 12 — > 24 and 18 — > 36. All the 
findings are similar to those obtained with s = 3/2.) 

For the s = 3/2 case with 3 lattice spacings and volume changes 12 — > 18, 16 — > 24 and 
24 — > 36 the results can be continuum extrapolated assuming a fit linear in a 2 /L 2 ; see the right 
panel of figure [I] and figure [2[ 

5. Conclusions 

In this work a running coupling scheme was investigated where the running scale /I is given by 
the linear size of the system [i = l/L. The main idea is to use the Yang-Mills gradient flow [|], ^, [|] 
but to adapt it to a finite volume setting [0]. In principle the original infinite volume construction 
can also be used to define a scheme where the running is via /i = 1 / a/87 where t is the flow time, 
however in this setup one would need to ensure that finite volume effects are fully under control. In 
contrast, incorporating finite volume dependence explicitly into the setup eliminates this problem. 
In addition, a small finite physical volume with appropriate boundary conditions for the fermions 
guarantees a gap in the spectrum even in the massless case. 

Our goal was to determine g 2 (sL) as a function of g 2 (L) at various lattice spacings. To that 
end we measured the renormalized coupling g 2 {L) and g 2 (sL) with the same set of bare couplings, 
m = and j8. In the next step we performed an extrapolation to the continuum limit, a 2 /L 2 — > 0, 
for the discrete j3 -function. 

For our main results we used simulations at 3 lattice spacings with s = 3/2 and c = 0.3. Renor- 
malized couplings between and 2.0 were studied and a controlled continuum extrapolation was 
carried out. For small renormalized coupling the universal continuum 1-loop result is reproduced 
and the MS 2-loop result is also consistent with our numerical results; see figure H This suggests 
that for our choice of s = 3/2 and c = 0.3 the coefficient a\ connecting the gradient flow scheme 
to the MS scheme g 2 = g 2 ^^ + a i( c )<?MS + • • •) i s probably small. 
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